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Abstract. We consider Cahn-Hilliard equations with external forcing terms. Energy decreasing and mass 
conservation might not hold. We show that level surfaces of the solutions of such generalized Cahn-Hilliard 
equations tend to the solutions of a moving boundary problem under the assumption that classical solutions 
of the latter exist. Our strategy is to construct approximate solutions of the generalized Cahn-Hilliard 
equation by the Hilbcrt expansion method used in kinetic theory and proposed for the standard Cahn- 
Hilliard equation, by Carlen, Carvalho and Orlandi, 1111 . The constructed approximate solutions allow to 
derive rigorously the sharp interface limit of the generalized Cahn-Hilliard equations. We then estimate the 
difference between the true solutions and the approximate solutions by spectral analysis, as in pQ. 
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1. Introduction 

In this paper, we apply an alternative method to matched asymptotic expansions, developed by Carlen, 
Carvalho and Orlandi, in [TT|, which allows the study of the sharp interface limit for the generalized Cahn- 
Hilliard equation, and derive higher order corrections to this limit. The method is based on the Hilbcrt 
expansion used in kinetic theory; we refer to [llj where the analogy is explained. We start by recalling some 
back ground regarding the Cahn-Hilliard equation and the results obtained in 

1.1. The Cahn-Hilliard Equation. Let f2 be a bounded domain in K 2 . The restriction of the analysis to 
two dimensions is made only for simplicity. Let m = m(x, t) be an integrable function on £7 which represents 
the value of a conserved "order parameter" at cc in SI at time t. The order parameter is conserved in the 
sense that J^ l m(x,t)dx is independent of t. Therefore, the evolution equation for m can be written in the 
form 

d t m(x, t) = V ■ J(x, t), 

where the current J is orthogonal to the outer normal of the boundary of f2. We take 

J{x,t) = a(m(x,t))V n(x,t), 

where a(m) is the mobility and is the chemical potential of x at time t. The mobility is positive and 

the chemical potential is defined as the L 2 (f2) Frechet derivative of a free energy functional T: 

K x ) = T- 0)- 
dm 

The simplest and most familiar example is the so called Cahn-Hilliard equation. It results by setting 
a(m) := 1, i.e. constant mobility, and 



F{m) :=\f |Vm(.*)| 2 d* + i f ( 



(m 2 (x) - l) 2 dx. 
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This leads to the evolution equation 

d t m(x, t) = A (-Am(x, t) + f(m(x, t))) , 

where 

(1.1) /(to) = to 3 - TO. 

Different choices of / can be made, provided they arc derivatives of a double well smooth enough potential 
with equal absolute minima. If m(x,t) is a solution of this equation, then 

A.F(m(-,i)) = - f \J(x,t)\ 2 dx, 

and thus, evolution decreases the free energy. The minimizers of the free energy are the constant functions 
to = ±1. These minimizers represent the "pure phases" of the system. However, unless the initial condi- 
tion too happens to satisfy J n uio(x)dx = ±|f2|, these "pure phases" cannot be reached due to the mass 
conservation law. Instead, what will eventually be produced is a region in which to rj +1 while m ~ — 1 
in its complement, with smooth transition across its boundary. This phenomenon is referred to as phase 
segregation, where the aforementioned boundary consists the interface between the two phases. If we "stand 
far enough back" from ft, all we can observe is the interface's shape since the structure across the interface 
is placed on an invisibly small scale. 

The evolution of m in time under the Cahn Hilliard equation, or another equation of this type, drives a 
very slow evolution of the interface. More specifically, let e be a small parameter, and introduce the new 
variables r and £ by 

t := e 3 t and £ := ex. 

Then of course it follows 

d t = e 3 d T and d x = ed^ . 

Hence, if m(x, t) is a solution of the Cahn-Hilliard equation and we define rn £ (£, r) := m(x(£), t(r)) then we 
obtain 

(1.2) 9 T TO e (C, r) = A e (- £ A c to e (£, r) + ^f(m £ (t, T ))j . 

If we think of e as representing the inverse of a large length scale, the variable £ will be dimcnsionless. The 
dimensionless variables are "slow" and the original variables "fast" for small e. In what follows, we keep 
the notation £ for the slow spatial variables, but we drop the use of r and replace it by t for convenience. 
One should just bear in mind that now we are looking at the evolution over a very long time scale when e is 
small. For the reasons indicated above, it is custom to consider initial data TOo(£) that is —1 in the region 
bounded by a smooth closed curve Tq in £7, and +1 outside this region. At later times t there will still be a 
fairly sharp interface between a region where m(£,t) ~ +1 and a region where m(£,t) ~ — 1, centered on a 
smooth curve T t . One might hope that for small values of e, all information about the evolution of m e {^,t) 
is contained in the evolution of the interface T t . This is indeed the case as shown in |llj . To explain the 
method used in [11], let Ai denote the set of all smooth simple closed curves in fi. As we will explain in 
Section 2, Ai can be viewed as a differcntiable manifold. A vector field V on Ai is a functional associating 
to each r in Ai a function in C°°(T). This function gives the normal velocity of a point on T, and thus 
describes a "flow" on Ai . We may formally write 

(1-3) ^r t = V(T t ). 

Now, given a flow on Ai, we can produce from it an evolution in C°°(Q) through the following device: Let 
to be any function from Ai to C°°(Q). We write to(£,T) to denote m(T) evaluated at £ G f2. We can then 
define a time dependent function m(£, £) on Q by 

(1.4) m(t,t):=m(t,T t ). 
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Notice that time dependence in m(£, t) enters only through the evolution of IV Now if, for small e and sharp 
interface initial data, all of the information about the evolution of solutions of the Cahn-Hilliard equation 
were contained in the motion of the interface, then one might hope to find a vector field V on A4 governing 
the evolution of the interface, and a function m from M. to C°°(tt) so that (|1.4[) defines the corresponding 
solution of the Cahn-Hilliard equation. 

In [11], a result of this type has been proved. More specifically, a sequence of vector fields Vq, Vi, V%, ■ ■ ■ 
defined on Ai was constructed such that the interface for the solution of (|1.2j) satisfies (|1.3|) for V := 
y"^. n e J Vj . It turned out that the leading term Vq is the vector field generating the Mullins Sekerka flow, 
as one could expect from the pioneering work of Pego [26] made rigorous by Alikakos, Bates and Chen [T|. 
In these papers the approximate solutions were constructed by using matched asymptotic expansions which 
give no information on the higher order corrections to the flow. The approach introduced in |llj unable to 
determine at any given order the velocity of the flow. 

Let us fix a number S > that will later be interpreted as a "surface tension', denote by K(£) = K(£, T) 
the curvature at £ £ T and by v the unit outward normal either to dfl or to T. Further, for each T in M, 
let \i be the solution of 

(1.5) A M (£)=o for £eO\r, 

subject to the boundary conditions 

(1.6) p(0 = S f JiT(0 - on T, d vt x = on 8Q, 

where |T| denotes the arc length of T and d v is the outward normal derivative to <9f2. Now define Vo(r) as 
the real valued function on T given by 

(1-7) K>(£,r):=i[cU] r (0 Cer, 

where the brackets on the right-hand side denote the jump of the normal derivative across V. In this way 
one defines a vector field on M. which generates a flow known as the Mullins-Sckcrka flow. For the local 
existence of a unique smooth solution of the free boundary problem (|1.5[) . (| 1 .€>[) and (|1.7[) see Chen,[T3] and 
Escher and Simonett, [T5]. As it is well known, the Mullins-Sekerka flow conserves the area enclosed by T t 
and decreases the arc length of T t - 

The higher order terms in ^j—Q^Vj are more complicated. In [TT], V\ which is the next correction to 
Vq was computed and a general technic of calculating all the higher order terms has been presented. The 
description of V\, like that of Vq, is in the context of potential theory. 

1.2. The e-dependent generalized Cahn-Hilliard Equation. We consider the generalized Cahn-Hilliard 
equation of the following type 

(1.8) d t m £ (g, t) = A ( - eAm £ (g, t) + t)} - G 2 ($; e)) + G 1 (^; e), £ in 0, t> 0, 

where A is the Neumann Laplacian operator on f2. The terms G\ and G 2 may depend on time also. In the 
present analysis, we shall consider the case where G\ and G 2 depend only on £ since our aim is to explain the 
main strategy in the simplest interesting setting. As it will be clear in the sequence, the proposed method 
is suitable for the time dependent case as well. 

The term G 2 in (|1.8[) models general external fields, see [23J [22] - In [22] the authors apply the Kawasaki 
exchange dynamics to derive a modified Cahn-Hilliard equation where G 2 describes the external gravity 
field. The free energy- independent term G\ may describe an external mass supply, cf. [22], or [3] where G\ 
was defined as a deterministic Gaussian function. Such a model is described for example in [3] , in order to 
model spinodal decomposition in the presence of a moving particle source, as a mechanism for the formation 
of Liesengang bands. In addition, G\ was introduced as a conservative white noise of thermal fluctuations 
cf. [53] or [TT] (Cahn-Hilliard-Cook model). Existence and uniqueness of solution for the stochastic problem 
was established in [151 1101 [B] while dynamics and stochastic stability were analyzed for the one-dimensional 
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case in [5]. Furthermore, the interface stochastic motion and singular perturbation has been studied for 
many related models like Allen-Cahn or Ginzburg Landau and phase- field models, cf. for example [20[ 14]. 
Integrating (|1.8[) over we get 

(1.9) d t (J^m e (£,t)d^) = J^8 V (- sAm E + - G 2 {s-,e))ds + J^G^e)^. 

Therefore, m e (£, i)d£ is not conserved unless the second member of (jl.9p is null. Generally, due to the 
presence of the external force field G2 and the external mass supply Gi, a free energy decreasing is not 
expected. For a mathematical analysis of the problem when G2 = and Gi is in L 2 (f2) cf. {27 ] fT6], 
An equivalent system formulation of (|1.8j) is the following 



(1.10) d t m e ^,t) = A^,t) + G 1 ^;e), 

(1.11) f(£,t) = -eAm^^t) + -/(m e (£,i)) - G 2 (£;e), 

£ 

where A is the Neumann Laplacian operator on fi. This representation will be used in our analysis. For the 
purposes of this paper we consider the e-dependent generalized Cahn-Hilliard equation (|1.8[) (and equivalently 
the system ([l.lOp . (|l.lljl ) supplemented with an initial condition 



(1.12) m e (£,0) = m^) 



1 on il Q 
1 on Oq, 



where £2q is the region of 12 enclosed by a smooth closed curve r and Oq = ^ \ (^0 u ^o) • Thus, we 
are assuming that the interface is already initial formed. Further we take the following Neumann boundary 
conditions 

(1.13) d v m £ = d v Am e = on <9Q. 

We assume that the forcing terms Gi and G2 are sufficiently smooth, and that d v G2 = on <9f2 so that 
(|1.13j) becomes 

(1.14) d v m e = d v [i e = d u Arrf = on dfl. 
We dot not require 

(1.15) / Gi(£eK = 0. 

J a 

Hence, mass conservation might not hold. The precise assumptions for the forcing terms G\ and Gi will 
be given in Section 2. For sufficiently smooth initial conditions and forcing terms G±, G2, there exists a 
unique classical solution of the generalized Cahn-Hilliard equation. The proof is analogous to that of the 
homogeneous case presented in [16]. 

Notice that if we write ()1.8|) in the original not scaled variables (x, t) the terms Gi and G2 are small 
perturbations of the standard Cahn-Hilliard equation. The term Gi in the original variables (x,t) is multi- 
plied by a factor e 3 and the term G2 by a factor e. The problem that we pose is the following. Take as in 
the homogeneous Cahn-Hilliard equation, initial data mo(£) like in (|1.12[) . Due to the presence of Gi and 
G2 the constant functions m £ = ±1 are not anymore stationary solutions of (|1.8|) . But we still expect that 
eventually at later times t there will appear a fairly sharp interface between the regions where m e (£, t) w +1 
and where m e (^,t) ss —1, centered on some smooth curve T t . We prove that this is indeed the case. We 
derive the motion of T t determining the vector field. It turns out that the leading term Vb in the vector field 
z^Vj (t\ N ^, governing the interfacial flow (see (|2.7jl ). is not the vector field generating the Mullins 
Sckcrka flow appearing in the sharp limit of the homogeneous Cahn-Hilliard equation. In fact, we obtain 

(1.16) Vo(;T t )^V {0) (;r t ) + (Vo)r t , 
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where 



/ v^ \ v ,T t )dS v ^o, 



and 

(1-17) (V )r t = WjJ ^ 6 "! *e[0,T]. 

Here, and in the following, we denote by dS v the element of the arc length along T or dd. We will indeed 
prove that, as e — > the singular limit of (|1.10[) and (|1 . llj) leads to the following moving boundary problem: 

Given a closed curve T° in that it is the boundary of an open set CIq C find a family jl^ £ M. : t £ [0, T] | 
and functions i) = fJ,(£,,T t ) for t £ [0,T] and £ G fl so that 
AM(f,t) = -Gi,o(0 Cen\r t , ie[0,T), 

/i(C,i)=25^(e,r t )-G 2)0 (0 on T t , ft, A i(.,t) = on dQ, *G[0,T), 
(1.18) 1 

^(■,r t ) = -[Mr t (0 ^ er *' *e(o,r), 
r = r°, 

where G\ o(£) := limGi(^;e) and G*2 o(£) : — limG2(£;£) and S > is the surface tension defined in (|4.32p . 

In [7] the authors applied formal asymptotics to analyze the sharp interface motion for generalized Cahn- 
Hilliard equations of the form (jl.8l) . The limit problem, which was formally derived in [7] , agrees exactly to 
(|1.18p which is rigorously proven in this paper. 
We immediately obtain for any t £ (0, T) 



(1.19) 



2/ V (r ) ,T t )dSr l = [ M r (v)dS v = - / A/ifo, t)d»7 = [ G^drj, 
Jr t Jr t Jn\r t Jn 



i.e (|1.17j) . Recalling that ^|O r J = J r Vb(?y, T^dS^, we obtain that the area enclosed by Tt is not conserved 
unless J n Gi t o(rj)dri = 0. Also, we have 

(1.20) ^|r t | = jT K-(77,r t )Vo(»7,r t )d5, = ^ Mv,r t )ds v + J^ v (v,r t )G 2 ,o(v)ds v ^ . 

Let us denote by /i ± (-,r t ) the restriction of /i(-,r t ) in f2^. It follows that 
2 / A t Vo(??,r t )d 1 S'^ = / fi + d v fi + dS, n - / yTd v n~dS v 

JY t JT t JT t 

(1.21) = - / div([iVfi)d£, - / <fw(//V/z)d£ 

Jq- Jn+ 

= - f |V M | 2 d£ - / ^A M d£ = - / |V/x| 2 d£ + f nGwffidZ. 
Jn JQ\v t Jn Jn\r t 

From these computations there is no reason to expect that ^|r t | is not positive. So, even in the case when 
the volume is conserved, i.e when J n Gi^o(r])dr] = the length of the curve does not decrease. The unknown 

|r t £ CI : t £ [0, T]| and ^ are coupled through the system (|1.18[) . However if the position and the 

regularity of the moving boundary |r t £ Cl : t £ [0, T]\ is known, the chemical potential fi is obtained by 
solving at each time t £ [0,T) the elliptic boundary value problem 
A M (£,t) = -Gx, (0 £eO\r t , t£[0,T), 

f i(U)=2SK(Z,T t )-G 2 ,o(0 on T t , d„(i(-,t) = on dSl, te[0,T). 
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In this sense we call a family {T(t);t G [0,T)} of surfaces a solution of (|1.18jl . To our knowledge there are 
no result regarding the existence and the uniqueness of solution for the moving boundary problem of the 
type (|1.18|) . A modified Mullins Sckcrka motion has been studied by [T7], but it differs from (|1.18[) either 
at the presence of the term — G\.q which is replaced in |17j by a specific function of time only, either at the 
presence of —G2.0 which does not appear in |17j . We think that a method similar to the one used in [17] 
might be useful to give existence and uniqueness of the classical solution of (|1.18[) . For the purposes of this 
paper, we assume that there exists a unique classical solution of the free boundary problem (|1.18l) . 

2. Notations and Main results 

2.1. Vector fields and flow on the curve space. Let Ai denote the set of all smooth simple closed curves 
in !lcR 2 . To discuss motion in Ai it is convenient to introduce local coordinates in the neighborhood of 
any given r G Ai. To this aim we define: 

Definition 2.1. Let K(^) = K(^,T) denote the curvature at a point £ G T for T G Ai. We define 

k(T) :=max\K($\. 

We denote by g?(£, T) the signed distance of £ G f2 from T. We define d < when £ is inside T and 
d > when £ is outside T. As long as T) < -j^j there is a uniquely determined point ijeT such that 
|r/ — £| = r); this is the point in T closest to £. Therefore, for any eo such that < Eq < /^py, let 

N(e a ) = N(e a ,T) := e : |d(£, T)| < e }. 

There is a natural set of coordinates in Af(eo). Given £ G AT(eo) we denote by p the diffeomorphism 
p : Af(eo) — > [— £o,£o] X T defined by = (d(£), s(£)) (whenever this does not cause ambiguity we omit to 
write the explicit dependence of N or d on T). We have that 

e = s(0+<M*(0), 

where v(s(£)) denotes the unit outward normal to T at s(^). For d E [— £o,£o] an d s G T let a(d, s) be the 
Jacobian of the local change of variables a(d, s) = dct^gjj—yl. A standard computation (cf. [21] appendix]) 

gives a(d, s) = YYi=i (1 — dKi(s)) , where Ki(s), i = 1, . . . , n are the principal curvatures at s G T, in the 
direction i. When n = 2 we have 

(2.1) a(d,s) = 1 -dK(s). 

In the sequel we identify functions of variable £ and functions of variable (d,s) in the domain Af(so). We 
denote by 

d 

z = — 

£ 

the stretched variable. 

The introduced coordinates in Af(eo,T) provide the means to give Ai the structure of a differcntiablc 
manifold and to study motions in this manifold, see [TTJ Section 2]. A vector field V on M. is a functional 
associating to each r in Ai a function in C°°(T). This function defines the normal velocity of a point on F 
and thus, describes a "flow" on Ai. More specifically, we may formally write 

(2.2) ^r t = V(T t ). 

We denote the lifetime T of the flow ()2.2[) , starting at T G Ai as 

(2.3) T = inf it > : k(T t ) < k \, 

where ko is any arbitrarily chosen positive number so that k(T) < kg < 00. If V(-, T) = K(-, T), the curvature 
at s G r, one obtains the curve shortening flow by curvature. When V(-,r) is given by (|1.7j) we have the 
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Mullins-Sekerka vector field, described in the introduction. When V(-,T) is given by ()1.16|) we have the flow 
characterizing the sharp interface motion studied in this paper. 

A given flow on M. produces an evolution in C°°(f2) through the following device: Let m be any function 
from M. to C°°(Ct); we write ra(£,r) to denote m(T) evaluated at £ € CI. Then a time dependent function 
m(^,t) may be defined on fl as follows: 

(2.4) m(£,t):=m(£,r t ). 

There is an obvious but useful decomposition of vector fields on M. . Given a vector field V on M. we may 
apply the decomposition 

(2.5) V(;T)=V^(;T) + (V) T , 
where 

(2-6) (V) T := ±-J V&T)dS 6 , 

and 

VM(;T) :=V(;T)-{V)r. 
Since (V)r is constant then by its definition is orthogonal to the constants in the L 2 (T) inner product 
i.e satisfies 

V (0) (r ? ,r)d5 r) =0, 



and therefore, it generates a volume preserving flow in the sense that for any t the area enclosed by T = T t 
is constant. 

Under the ansatz given below, in this paper, we derive separate equations for the components V^°)(-,r) 
and (V)r for each of the vector fields Vj. 

Ansatz 2.2. Let Vo, V\, V%, ■ ■ ■ be a sequence of vector fields on M. and m , mi, m^, ■ ■ ■ functions from M. 
to C°°(f2). For any given initial interface Tq in Ai and all N > 0, let t[ N ^ be the solution of 



(2.7) dit 



N-l 



(if) with rw=r . 



dt 

We define the function m^ N \^,i) by 

Ait T l ( 7V )^ N 

(2.8) mW(f,t) =m„(^^ + ]>>m,(£, rf>), 

£ 3=1 

and notice that m^ N ^(^,t) depends on t only through 
We set 

(2.9) m {z) := r(^-z)m(z) + (l - r (i-z))sgn(z), 

where m(z) := tanh(z/ v / 2)0 defined for any z € K is the unique solution of the Euler- Lagrange equation 

-m"(z) + f{m(z)) = 0, z e R, lim m(z) = ±1, 

z— »±oo 

and r is a smooth even unimodal cut-off function, r{u) = 1 for \u\ < ^ and r(u) = for u > 1. 
In addition, let 

a(c t( n ")\ 

(2.10) m^£ ) rW):=/ 1 ,(^^ lS (c J rW))+^(c,rW), e e n, ; = i,-,iv, 

^The explicit form of the solution is never used. We will use only its qualitative properties. 
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where hj are C°°(Q) functions equal to in fl \j\f(eo) and when d(^,V[ ) = 0. The functions (f>j, j = 
1, • • ■ , N are in C°°(f2), satisfy the Neumann boundary conditions on dQ and admit a global Lipschitz bound, 
independent of e, i.e. 

\\<t>j\\u P (Q) < C j = l,---,N, 
where C is a constant independent of e. 

Notational convention Below we denote by mW(^t) := mW(£, y[ N) ) and ^ N - l \^,t) := ^ N ~^ (£, T { t N) ). 
If there is no ambiguity we write T or T t for T^ . Trough what follows, we write C to designate a generic 
positive constant independent on s. Its actual numerical value may change from one occurrence to the next. 

Remark 2.3. The Ansatz \2.2\ must be modified when G\ and G 2 depend on time. The (|2.8[) should be 
replaced by 

jit r( N h N 
mW ( fl t) = to + £ tt, t, if )- 

3=1 

Notice that m^ N ^(^,t) depends now on t not only through t[ N . One can verify that the first order function 
to keeps to depend on t only trough T^' . 

2.2. Main results. We start constructing a function m( N >(£,Tt), for £ 6 £1 and £ <G [0,T] where T is the 
lifetime of (|2.7[) and show that it is an approximate solution of (|1.8[) . We make the following assumptions 
on the forcing terms G\ and G2. 

Al: Assumptions for G\ and G2. For any TV > 1 we require 

N-l 

Gifae) = E £ J G hJ (0 + e N G hN (Z,e), \G t , N (te)\ < C, 

(2.11) J =0 

G,., ■ C°°(fi) i = 1, - - - ,N-1, i = l,2, 

d v G 2tj =0 on 90, j = l,---,N. 

Remark 2.4. M^e require Gij € C°°(Cl) for j = 1, • • • , TV — 1 and /or i = 1,2 to avoid regularity problems, 
but this assumption can be relaxed. 



Theorem 2.5. Let N > 1 and G\ and G2 be as in assumptions Al. There exist vector fields Vj, j = 
0, • • • , (iV— 1) and functions mj , j = 0, •• • , N as prescribed in the Ansatz \2.2\ havina the following properties: 
Let T denote the lifetime of the solution of (|2.7[) in Ai . Then there is a constant Cm so that for all t < T 



N-l \ N-l 



(2.12) cW^tt,*) = A -eAm^tt,*) + - V e>G 2 A0 + E e j G ltj (0+AR^^,t) 



3=0 I 3=0 



< Ov£ 



where 

(2.13) sup R {N) {t,t) : 

£en,te[o,T] 

Finally, the sequences of vector fields and functions are essentially uniquely determined since given Vj for 
j < k then Vk is determined up to 0(e k+1 ), and similarly, given rrij for j < k then is determined up to 
0(e k+1 ). 

Remark 2.6. Ln Theorem \2.5\ and in the following, the symbol 0(s m ) denotes terms which are of order e m 
uniformly in all their variables. The qualified nature of uniqueness stated in this theorem is an indication 
that there will be choices to be made at every stage of the approximation. 
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The proof of Theorem 12.51 follows the main lines of the scheme introduced in and it is proven in 
Section 5. There, the complete result relating the solution of (|2.12[) and its sharp interface limit is given. The 
construction behind the proof is patterned on the Hilbert expansion of kinetic theory. We refer the interested 
reader to Subsection 3.2] where this connection is discussed. We first construct an approximate solution 
up to order N of the chemical potential /i e (cf. (jl.lOjO assuming that the left hand side of (|1.10[) is known 
and it is given by the Ansatz 12.21 This is done in Section 3. We, then, insert the constructed approximate 
chemical potential into . The approximate solution m W is determined provided certain compatibility 

conditions are verified. This is done in Section 4. Finally in Section 5 we construct (rh^ N \ /i^ - 1 )) where 
rh^ N ^ is an e N modification of mS N ^ and jl^^ 1 ^ is an e N ~ 1 modification of fi^ N ~^ and we show Theorem 



Let Vj, j = 1, • • • , N — 1, be the sequence of vector fields introduced in the Ansatz 12.21 According to 
(|2.5j) . wc split them as 

v j (-,T) = v / (0) (.,r) + (y J ) r , j = i,---,N-i. 

The term Vj is determined in Theorem 14.11 by applying the Dirichelet- Neumann operator in the context 
of potential theory, while the term (V^r, which is constant on V, is determined in Theorem 13. II 

As already explained in the introduction, the leading term Vb in the vector field "^2^=0 ^Vj y^t*^ 
governing the interfacial flow (cf. (12.71) ) is given by 

V (;T t )^V^(;T t ) + (Vo)r t . 

The family of curves ^Tt G M.,t g [0,T]| driven by Vb is the solution of the moving boundary problem 

(|1.18|) . For a curve r £ M. . the term Vq°^ (T) is determined, see Lemma 14.51 as a real valued function on T 
given by 

(2.14) K, (0) (e,r)-[a,MoAo] r (C) eer, 

where the brackets on the right denote the jump in the normal derivative across T. For each r € Ai, ^o,o,o 
is the solution of 

(2.15) A/*(£) = for £ 6 Q \ T, 
subject to the boundary conditions 

f T K(QdS 6 \ 1 / . e , / r So,o,o(Od5e 



(2.16) M (0 = S \K® - Ji |r y ' )+- (A,,o,o(0 - Ji ""gp " J on I\ d vl t = on dCl. 
The term B ,o,o ( cf - 6330 is g iven by 

(2.17) S ,o,o(0 = -2^,0,0(0 + ^2,0(0] 
where (cf. (081 1 



(2-18) Ao,o,o(0 = (V&> r 2 / G(£,r))dS v - / G^,r,)G h0 (r,)dr,, 

for G(£, rj) the Green function in f2, with Neumann boundary condition on dfl, satisfying the equation 
(2-19) =5{£-ti)- p. , 

so that 

(2.20) / G(Z,ri)dri= f G(tv)d£ = 0. 

Jn Jn 

The second step is to show that there is an actual solution of (|1.8[) close to the constructed approximate 
solution m^(-, •) whenever both of them start from the initial datum mg, see (|1.12p . This step for the 
standard Cahn-Hilliard equation (i.e. for G\ = G2 = 0) has been proven in the work of Alikakos, Bates and 
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Chen, PP, by application of spectral estimates. We use the spectral estimates as in pp. Namely the liner 
operator 

Cw = A (eAw - V(mW(t))? 

that one obtains linearizing p. 81) at m^ N \t), t <= [0,T]; the solution constructed in Theorem l2.51 is the same 
as in pp. The approximate solution has the requirements needed to apply the spectral estimates proven 
by [2] in two space dimensions and by [12] in arbitrary space dimensions and in more general setting. The 
assumptions imposed on the forcing terms G\ and G2 together with the assumption (|2.21[) are enough to 
have these terms under control. We state the theorem and we outline in the appendix the proof. 
Let p > and |j • be the usual norm in L p (Vl), then for T > we define the norm 

I f T \ 1/p 

\\u\\p,Q T ■= [ j o IMIp,I2*J ■ 

Theorem 2.7. Take N > 1, G\ and G2 as in Assumption Al. Further assume that 
(2.21) / G liJV (&e)de = 0, Ve>0. 



Letm^ N \t) forte [0,T] wh ere T is the lifetime of (|2.7I) , be the solution of (|2.12[) . Let m 6 be the solution 
of the generalized Cahn-Hilliard equation (|1.8[) supplemented by the boundary conditions (|1.14l) and having 
initial datum m e (£,0) = m'^^, 0), £ € fl. Then, there exists So > so that for all e <G (0,£o]; f or an U P a * r 
A > N > ^A±5, it holds that 

(2.22) ||m £ -mW|| 3 ,n T <£ A 

Remark 2.8. The result of the previous theorem coincides with the analogous result in L p (Qt) norm of 
Theorem 2.1 of [T] for the case G\ = G2 = in dimensions n = 2, since 3 = p — 2^g and since 

A>f = (n + 2)Bi±M 

Remark 2.9. Set Gi.n-i = G^jy— 1 + zGi.n and therefore G^jv = 0. Determine the velocity field V/v-i 
replacing G^/y-i with Gi.n-i- In this way the condition (|2.21[) of Theorem\2.7\ is trivially satisfied. 



3. The Construction of the Approximate Chemical potential 
In this section, we apply potential theory to show the following result. 

Theorem 3.1. Take N > 1 and G± as in Assumption Al. Let r[ , t <G [0,T], be the solution of (|2.7[) 
in Ai, with T its lifetime, see (1231) . Let m^ N \;T\ N) ) be as in the Ansatz \2.2\ There is a unique way 
to determine the (Vj)(T^), j = 0, ■ ■ ■ ,(N— 1), such that there exists a unique (up to a constant in 
expansion 

N-l 

(3.1) ti {N - 1] (t t) = e >^> f ) m n x (°> T )' 

i=0 

with 

N-i 

(3.2) = A/i^-^^.t) + XI £^^-(0 + i?i(C,i,e) (0,T), 

with Ri given in (j3~TIj) . Fwrt/ier At (JV_1) (-, £), /or f G (0,T), is a C°°(f2) function satisfying the Neumann 
homogeneous boundary conditions on dfl, 

(3.3) sup \Ri(^,t,e)\ < C(T)e N ~ 1 , 

5,tenx[o,T] 
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and 
(3.4) 



sup | f R 1 (^t 7 e)d^\<C(T)e N , 
te[o,r] Jn 



where C(T) is a constant independent of e. Moreover, the terms fit appearing in <|3 . 1 1) are specified by (|3.18|) . 
([3^27)1 and ([338)) below. 



We look for a function ^ N ^ from M to C°°(n) having the form 

N-l 

(3.5) 



^- ^r) = J2 e*Mi(^r) C e n, &w - o on an, 



,(JV- 



i=0 



where /i,;, i = 0, • • • , N — 1 are functions to be determined. We insert into (|1.10[) the function mS N \ given by 
the Ansatz |2~21 and ^ N ~^ given by (|3.5p , where both are evaluated at Y := y\ N \ for y[ N ^ the solution of 
(|2.7[) . Therefore, we obtain (AT — 1) Laplace equations for fj,i(-,Y[ N '), i = 1, • • • , (iV — 1). The compatibility 
conditions are needed in order to solve these equations and determine (Vj)(Y\ N ') for j = 0, • • • , (AT — 1). 

When differentiating m^ N \-,Y\ N ') with respect to t we need to take into account that m^ N ^ depends on 
r t through a fast and slow scale. The fast scale brings a factor e -1 . 

Definition 3.2. Let m be a function from M. to C°°(tt) of the type h(^ d ^ r ^ , s(£, Y)j and V be a vector 
field on M. . We define 



(3.6) 



D v m(£,Y) := -h! 



i t ,M£,r) 



,s(Z,Y))V(s(0), 



where the prime indicates the derivative of h with respect to the first variable z = d ^ ,r ) . 



In addition, for any Wn '■= X^Lo* w ith ^bi ' ' ' j Viv— l vector fields on M., we define 



N-l 



(3.7) 



D WN m(^Y) := £ e j D Vj m(Z,Y). 

3=0 



Note that by the orthogonality of V^d with respect to the surface there is no contribution in (|3.6|) from 
s(£, r). Therefore, cf. [TT] for the detailed computations, differentiating (|2.8|) with respect to t, applying the 
chain rule at the right-hand side (here the velocity will appear since rrij are defined on Y) and then using 
flUEl , ([376]1 and ([37f]> . we arrive at 



d t m^ =D Wn (m (jv:) ) = D Vo m + e [D Vl m + D Vo mi] 



(3.8) 

where 
(3.9) 



e 2 [Dv^mi + D Vo m 2 + D V2 m ] 



JV-1 



N-l 



i=0 



7? 



2V-1 



i=0 



+ 0{e N ). 



The mo is the function defined in (|2 .9[) and the term _E = £7(£, t, e) is obtained by differentiating r( ■ ) , 

. £ ° 

the unimodal function appearing in the definition of mo, with respect to the velocity field. E is given by 

[N-l "I 

XVVi(a(0,t) {m-[l 



(3.10) 



t, = — r 



to 



£0 



M£,r t )>o} - I {d(«,r t )<o} 
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Note that E is exponentially small since r' is different from zero only for || < \z\ < while fh converges 
exponentially fast to ±1 as z — > ±00, jTT]. Taking into account (|3.8I) and (|1.10j) we obtain a set of N 
equations for the /ij, i = 0, • • ■ , N — 1. 

Zero order term in e: 



(3.11) 



D Vo m = ivbm' = A/ao + Gi,o 



9 y /Lto = on 90, 
where G± is the zero order term in e of G\ . 



for 



(efi, t€[0,T}, 



First order term in e: 



(3.12) 



D Vl rriQ + Dv n mi = Afi\ + G\.i 
^Mi = 011 <9£1. 



for 



ie[0,T], 



n-th order term in e (n < N — 1): 



(3.13) 



— A/x„ + Gi n 



d v Hn = on dfl. 



for 



(efi, ie[0,T], 



Remainder term: The remainder term, see (|3.9p and (|3.10[) . is given by 
(3.14) fl x (f , t, e) = e^Gi.jvCf , e) + i?jv(£, t, e) + i, e). 

Since the derivative in Rn (cf. (|3.9p ) brings down a factor e -1 then this yields easily the next estimate for 

Ri 

(3.15) 



sup |i2i(e,*)| < ccr)^ -1 - 

(5,t)Gf2x[0,T] 



Further, one gains an extra power of e when integrating since the terms of order e N 1 have support in 
AA(e ) 

(3.16) sup / IR^trfldZKCffle". 

te[o,T] Jo. 

In the sequel, we prove existence and uniqueness (up to a constant) of solutions of the equations obtained 
so far at different orders. In Lemma 13.31 and in Lemma 13.51 we consider the first and second order equation 
respectively. Finally, in Lemma 13.61 we outline the proof for solving the equation to a generic order. In the 
next lemma we write in an explicit way the dependence of the mean velocity on s. This is done in order to 
get easily the leading velocity field governing the interfacial flows, see (|2.7p . 

Lemma 3.3. Under the conditions 

V (;T t ) = V {0 \;T t ) + (Vb) rt [l + co(e)], 



(3.17) 



V^faTjdSn^O, (V )r t 



1 

2|T 



t\ Jn 



G lfi (v)dv t€[0,T], 
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where co(e) defined in (|3 . 24[) goes to zero exponentially fast as e — > 0, there exists a unique solution (up to 
constants in of p. lip given by 

(3.18) r t ) - /i 0>0 (e, r t ) + Ao(e, r t ), 

where 

(3.19) Mo,o«,r t ) = / G(t, V ) (-m' f^hM) V^°\a(r,),T t )) d v + c (t). 



Here, co(t) is a constant (in to be determined, and 

(3.20) mo(£, r t ) = (V )r t [1 + co(e)] ^ G(£, »?) QmJ, (^^) ) dn - J G(£, v)Gi,o(vW 
The term no is in C°°(f2) /or any i € [0,T]. 

Proof. Because (9„/io = on 90, the solvability of (|3.1ip requires that for all t E [0,T] 

(3.21) J f^m' f^^) V (s( V ), t)\ dn J Gx^dr, = 0. 

In two dimensions (note that in three dimensions there will be extra terms), by using local coordinates it 
follows that 

\< {^P) Vo(s(v),t)) IV = l i6o) (^) V (s(r,), t )) d V 

m' (z)Vo(s,t)e(l — ezK(s))dsdz 

£ Jr J—Sl 



1 

(3.22) 



m' (z)Vo(s, t)dsdz — e / / zm' (z)K(s)Vo(s, t)dsdz 

m' (z)V (s,t)d S dz = 2(l-e-^) f V (ji,T t )dS n . 

The last line holds true since m' is even and exponentially decreasing. Replacing now (|3.22[) in (|3.21l) . we 
obtain 



(3.23) G l>0 {ri)dr, = 2{l-e--) V (r),T t )dS v . 

Jn JT t 

Taking into account the splitting Vb = V^ ' + (Vb)r[l + co(e)], (cf. the first and second equality in (I3.17[) ). 
and using (|3.23p . we arrive at 

V (ri,r t )dS v = \r t \(V )r[l + c (e)] = hl + c (e)) f G lfi {rj)dr}, 

where 

(3-24) co(e) := 6 * . 

(1-e O 

This forces to take Vb satisfying the third relation of (|3.17p . By potential theory, once the compatibility 
condition is satisfied, the solution is given by ()3.18j) . □ 

Remark 3.4. Note that (Vb)r t an d Ao(£j-Tt) are completely determined once we know m®. The quantity 
Mo.o^rt) depends on Cq{£) and . These quantities will be determined when proving Theorem \4-l\ 
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3.1. The first order term in e. For the derivation of the first order correction we need to prove the 
solvability of equation (|3.12l) . 

Lemma 3.5. There exists a unique (up to constants in solution [i\ of (|3.12[) provided that 

(3.25) ^(r^y/'VtHWr*, 

with 
(3.26) 



/ ^ 0) (77,r t )dS, = 0, Vie[0,T], 



and (Vi)r t chosen according to 

{Vi)v t 



■[l + c (e)] 



/ Gi,i(Ode-&i(*) 



where b±(t) is defined in (|3.30[) and cq(e) in p.24[) . The solution is given by 

(3.27) Mi(£,t)=Mi,ote,t)+£i&t). 
where fix is defined in (|3.35|) . while 



(3.28) 



VH, (£,t)= I ( ~< 



1„, r^lV^K^^ld. + cUO- 



Here, c\{t) is a constant (in to be determined. In addition, the solution is a C°°(f2) function for any 

te [o,T]. 

Proof. The solvability of (|3.12[) requires that 

(3.29) / [^TOo + ^milde- / Gi,i(Ode = 0, 

for any t G [0, T]. Here, we are assuming that mi, mo and Vo are already determined and so we define 

(3.30) h(t) = / D Vo mid£. 

Jn 

Proceeding as in (|3.22|) we obtain 



(3.31) 



D Vl m d£ = 2(1 -e~) V 1 (r,)dS r ,. 



Taking into account the decomposition p.25[) and relation p.26[) we have 

/ V 1 (v,Tt)dS v = \Tt\(V 1 )r t . 

Therefore, relation p.29[) is satisfied if 



(3.32) 



0"i 



i r. 



2T t 



■[l + co(e)] 



/ Gi,iKK-&i(*) 
Jn 



where co(e) is defined in (|3.24|) . This determines (Vi)r t) the projection of Vi(T t ) onto the constants. The 
solution of (|3.12|) exists and it is given by 



(3.33) Mi(f.*)= / Gfarj) D Vl m + D Vomi dry- / r))G ltl (r))dr) + c x (f). 

Since we shall use the decomposition (|3.25|) . it is convenient to write 
(3-34) Pi(&*)=^i,o(£,t) 
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where ni,o(£,t) is given in (|3.28j) and fii is defined as follows 

(3.35) , t) = J G(£, ^Dv.rmdq + (V x )r ( j G(£, r,) Qm^ (^lEil 



d77- / G^G^dr). 
Jn 



□ 



Lemma 3.6. The solution fij(-,t) of (|3.13[) for 2 < j < N — 1 exists and is unique (up to constant in £,) 
provided that 



(3.36) 



(3.37) 



^.(r t ) = ^°)(rt) + <^>r sl 
J v}°\s,T t )ds = o, v<e[o,r], 



and (Vj)r t is chosen according to (|3.44[) . It is given by 
(3.38) 
where 



(3.39) 
and 

(3.40) 



V^ >( a (f?),*))d^ + Cj -(*), 



( 7™o ^ 



= / G(e^) 



3-1 
n=0 



(1// 



+ (^)r t / G(Z,r,) (-m' 



1 , / d(ri,T t 



The solution fij{-,t) for t £ (0,T] is a C°°(f2) function. 

Proof. The proof is analogous to the proof of Lemma 13.51 The solution exists if 



(3.41) 


/ 


j 

^2 Dv n mj- n 






Jn 


n=0 


Jn 



for any t £ [0, T]. Here, Dv n mj- n for n = 0, • • • , j — 1 arc determined and so, we define 

/" P'" 1 

(3.42) 6 3 -(t)= / VlVnTn,--, 

[ n=0 

Requiring (|3.36[) and (|3.37[) we obtain 



(3.43) 



D v .m de = 2|r t |[l-e--](^)r t . 



Hence, to fulfill relation p.41j) . we must take 

1 



(3.44) 



(Vj)r t 



[l + c (£)] 



G^iOdH-bjit) 



3Ut 2\T t \ 

where co(e) is defined in (|3.24p . This determines (Vj)r t , the projection of Vj(T t ) onto the constants. It still 
remains to determine the orthogonal part Vj°K The solution of (|3 . 1 3[) exists and is given by p.38[) . □ 

3.2. Proof of Theorem 13.11 From Lemma T3.3[ Lemma T3. 5 1 and Lemma T3. 61 we have that /x^" 1 ) satisfies 
by construction ()3.2[) . The remainder R\ is defined in f|3 . 14[) and estimated in (|3.15[) and (|3.16f) . The term 
fj,^^ 1 ' (• , t) for t £ [0, T] satisfies the homogeneous Neumann boundary conditions by construction. Thus, 
Theorem 13.11 holds true. □ 
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4. Derivation of the equations for 

Next theorem assures the existence and (essential) uniqueness of the functions rrij, j = 0, ■ • • ,N, having 
the properties required in the Ansatz 12.21 Existence and uniqueness are obtained provided that a certain 
compatibility condition is satisfied. This determines vj°\ the orthogonal part of the velocity field. 

Theorem 4.1. Take N > 1 and G 2 as in Assumption Al. Let T be the lifetime of the solution of (|2.7|) 
in M. Let *)(•, t), t G [0,T] be the function constructed in Theorem \3.1\ Then it is possible to choose 
the vector fields Vj so that there exist rrij, j = 0, • • • ,N, having the properties prescribed in the Ansatz UTB 
such that 

N-l 

(4.1) ^- 1 )(e i t)=-£AmWfet) + -/(roW(C,t))- £ e*G 3l3 -(0 + ik(£,t,e) in Qx(0,T], 

with i?2 given by (|4~7T|) . Further, mW(-,t) forte [0,T] is a C°°(fi) function that satisfies the homogeneous 
Neumann boundary conditions and 

(4.2) sup sup \R 2 {£,t,e)\ < Ce N . 

fefi te[o,T] 

Finally, the choice of the term vj- * 1 is specified by the equations (|4.3ip . (|4.5ip and (|4.63[) given below. 

In p. lip we insert at the left-hand side the already determined function /i^ -1 ^-, ) and we obtain 

(4.3) m (JV-1) (£ *) - -sAmfo t) + * )} - G 2 (£; e) in X (0,T). 

e 

Then fx^ N ~^ is written in terms of m^ N \ chosen according to the ansatz. Here, we prove that there exists 
a unique way to find the function m^ N \ having indeed the property required in the ansatz and satisfying 
equation (|4.3[) in a certain sense (to be specified in the sequel). 

The existence of any rrij, j = 0, • • • , N is obtained provided that a compatibility condition is satisfied. 
This compatibility condition forces us to define properly vj°\ j = 0, ■ ■ ■ , (N — 1). 

We distinguish two main steps: 

• Step 1: Determination at any order of the equations. This is carried out in the Subsection \4-l\ 

• Step 2: Analysis of the equations derived at step 1. This will be done in the Subsection \4-2\ 

In the present and in the next subsection T t is kept fixed, therefore, for the sake of a simpler notation we 
drop the subscript t, except of the cases that a subscript use may add some clarity in our arguments. 

4.1. Determination of the equations for rrij, j = 0, • • • N. To separate the fast and slow scale of rrS N ^ 
near the surface T, we write the Laplacian in the system of local coordinates introduced in Subsection 2.1. 
The expansion in e of the Laplacian written in this coordinate system is reported in the Appendix, subsection 
A. 2. We then match the right and left terms of the equations having the same power of e, distinguishing the 
case where £ € A/"(eo) from the one with £ £ tl\Af(eo). We therefore, get at any order two sets of equations: 
one for £ £ Af(eo) an d the other for £ G £1 \ J\f(eo). Since the interface separates ft in two regions we will 
distinguish those £ G Q \ Af(eo) which are inside T from those £ G £1 \ Af(so) which are outside T. 

Taking into account formula ([6. lip in the Appendix, denoting by ' the derivative with respect to z, and 
by a n , b n , c n the quantities defined in (|6.12p , after simple, however lengthy computations we obtain the 
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following identity 



JV 



: 2 AmW(z, s) = I fh"(z) + £ e n K(z, s) + a n {z, s)m') 



n=l 
JV n-l 



JV 



(4.4) 



^n=2 i=l 
JV n-3 



n=3 



(z, S )^(z, S )|>+£ 2 A 



n— 4 2=1 



^ b n -i(z, s) — h l {z, s) 

;=i 

TV 



with 
(4.5) 

(4.6) 



sup L4(£,t, E )| <C(T), 
(f,t)enx[o,T] 

sup / de|A(e,*,e)| <eC(T), 

te[0,T] Jf2 



(4.7) Ek&e) = e 2 Ar(^Il) |m(^|^) - [l MC , r) >o} - I {d(5 ,r)<o}] 



m( ) 



and 
(4.8) 

t 

Let us now define fi such that 
(4.9) /(mW) = /(m ) + / / (m ) 



lim sup |S 1 (f,t,e)|=0. 
e_K) K,()enx[o,T] 



' JV 

E 



£ TOj 



^£7i( m 0i m i.",m«-i) + £ Ar+1 S/v+i (■,£), 



(4.10) 



sup \B N+1 %t,e)\<C. 
fen,te[o,T] 



One can easily obtain /,; for any i = 2, • • • , JV by using Taylor expansions up to JV — 1 order for / around 
mo and collecting then the terms of the same power of s. 

We insert (|4.4[) and (|4.9[) into (|4.3|) and equate terms having the same order in e (when estimated with 
the L°°(f2) norm) obtaining, this way, two equations at any order, one for £ £ fi \M(eo), and the other for 
£ € Af(eo)- The equation for £ £ O \ jV(£o) will determine fa which are the slowly varying terms, while 
the one for £ £ Af(eo) will determine hi i.e. the rapidly decaying terms. When deriving the equations for 
£ £ M(sq) then terms of the type Afa(sz, s) appear. The fa are C°° functions, since they are proportional to 
/ij, and they have the same type of singularity in £ when differentiated in £. So, the terms £ ,l+1 A</3„_i(£z, s) 
are of order 0(£ ra ), and therefore, we write them in the e n order equation. 

In the following we will use the notation /(±1) (or /(m(±oo))) in f2\ Af(eo). This refers to the fact that 
r separates Q in two sets, i.e. fi + which is the part outside T and the part inside T. Therefore, ± refers 
to different regions of (f2\A/'(£o))nJ7 ± . It is convenient to write the external potential G^i, for i = 0, • • • , JV 
in local coordinates when £ £ jV(£o)- We therefore identify G2,i(£) = G2,i(d(£, T), s(£, T)), for i = 0, • • • , JV. 



Zero order term in e: Matching gives 

^£o 



(4.11) 







r( — z )m (z) + f(m (z)) for z€ 

V£n / 



£o £o 

£ £ 
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and 

(4.12) /(m (±oo)) = for £ e fi\JV(eo). 

First order term in e: Again matching gives for z £ [— ^r, ^r] and s £ T 

(4.13) Mo(ez,s) =-[til(z,s)-K(s)m' (z)} + f'(mo)[hi(z,s) + <i> 1 (ez,s))-G2fi(ez,s) 
and 

(4.14) Mo (0 = /'(m(±oo))^(0-G 2 , (£) for ^£n\N(e ). 

Second order term in e: 

(jLi(ez,s) = - [h'2(z,s) - K 2 (s)zm' (z) - K(s)h' 1 (z, s)] 

(4.15) 

+ / (m (z)) [h 2 (z,s) + 4> 2 {ez,s)] - eA0i(ez,s) + f 2 (m ,m 1 )(ez, s) - G 2 ,i(ez,s), 

(4.16) m (O = /'(m(±oo))02(O + /2(m(±oo)>i(C))-eA0 1 (O-G 2 ,i(O for £efi\^(e ). 
More explicitly the f 2 term for £ £ M(sq) is given by 

f 2 (m ,mi)(sz,s) = -f"(m (z)) [h\(z, s) + (j>l(ez, s) + 2<f> 1 (ez, s)h 1 (z, sj] , 



and by 

for£e n\Af(e ). 



/ 2 (m(±oo),fc(0) = i/"(m(±oo))^(0, 



n-th order term in e (3 < n < N): 

n-l 

Hn-iisz^) =- h'^(z,s) +a n (z,s)m' (z) +^2a„- i (z,s)h' i (z,s) 



n-2 ^ 2 «- 3 ^ 

(4.17) + ^ b n _i(z, s)-^hi(z, s) + 1{„> 4 } ^ c n-i(z, s)—hi(z, s) 

i—1 i—1 

- eA<j> n -i(ez, s) + f'(m ) [h n (z, s) + <p n (ez,s)] 

+ f n (m ,m 1 ,m 2 , .., m n -i){ez, s) - G 2 , n -i(ez, s), £ £ Af(e ), 

(4.18) M n-i(0 = -eA^-itf) + /'(±1)^„(0 + /n(±l, 01, fc, ... 0n-i)(O - G 2 ,„-i(0, 
for£e O\AA(e )- 

Remainder term: The remainder R 2 (£,t,e) = R 2 is the following: 

(4.19) £j R 2 = e" +1 G 2>JV + e^+M + e N+2 /\<p N +E l+ e N+1 B N+1 . 
From (|4"3|) . (f^Tgj) and (|4~TU)) we obtain 

(4.20) sup |j2a(f,t,e)| < Cfe*. 

(C,t)enx[o,T] 

Remark 4.2. Since, for i = 0, • • • , -/V, we required d l/ G 2 j = on dfl then the fa constructed in Theorem 
\3.1\ satisfy d v fa = on dVL; to obtain d v rrii = on <9f2 it is enough to have d u 4>i = on dfl. 
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4.2. Analysis of compatibility conditions. Our aim is to analyze the equations obtained so far in the 
previous subsection. The strategy is to find first at each order in e the slowly varying part <pi by solving 
the equations for £ G f2\ N{eq, Y\ N ' ). Then, we extend globally in fl and determine the rapidly decaying 
part hi by solving the equations in £ € jV(4p, T\ N '). However here, in order to continue to arbitrary order, 
it is convenient to modify the way we extract the compatibility condition required for solving the equation 
for hi. The modification is to add and subtract to each order a term of lower order e z+1 ai(s, T)fh'(z), 
with cti(-,T) <G C°°(r). Adding and subtracting terms does not change, of course, the total quantity but it 
modifies the equation we obtain at each single order. In the following, for the sake of short notation, we 
write ai(s) = ai(s,T). 

At each order i > 1 in e the associate function m, splits into two parts. The first is the function <pi defined 
globally in Q and satisfying Neumann condition on the boundary of Q while the other part is the function 
hi which differs from zero only in a tubular neighborhood of T, Af(eo,^), and is exponentially decaying to 
zero far from T. 

The zero order term is different in the sense that mo far from the interfaces relaxes exponentially fast to 
±1. We first state the following lemma, which is taken from [1]. We shall use this lemma to determine the 
condition for solvability of equations of the type (|4.23p , where £ is the operator on L 2 (R) defined by 

(4.21) Cg{z) = -g"(z) + f'{m{z))g{z). 

The operator C is self adjoint on L 2 (M.) and has a null space spanned by in'. Therefore, the condition for 
solvability of Chi = g is 

(4.22) / g(z)m'(z)dz = 0. 

Lemma 4.3. [see [1]/ Let A(z,s,t) defined for z £ 1, s 6 T, ( £ [O,? 1 ]. Assume that there exists A ± (s,t) 
such that A(z,s,t) — A ± (s,t) = C(e~ Q l z l) as \z\ — > oo for seT and t € [0,T]. Then for each s € V and 
t e [0, T] 

(Cw)(z,s,t) = A(z,s,t) for z e K, 
< ~ 4 ' 23 ' ) w(0,s,t)=0, w(-,s,t) G L°°(R), 

has a solution if and only if 

(4.24) A(z,s,t)m'(z)dz = for all s e T, ie[0,T]. 

Jr 

In addition if the solution exists, then it is unique and satisfies 



(4.25) Di 



, A ± (s,t) 



C(e" a|z| ) as \z\^oo, for ^ = 0,1,2. 



Furthermore, if A(z, s,t) satisfies 

D" s l Dl l D e z [A(z,s,t)~A ± (s,t)] =G(e- Q|z| ), 

then 

DTD-D e z L(z,s,t) + ^^=0(e-<*^), 

for all m — 0, 1, ■ ■ • , M , n = 0, 1, • • ■ , N , and I = 0, 1, • • ■ , L + 2. Further, since C is a preserving parity 
operator, the solution w(z,s,t) is odd (even) with respect to z if A(z,s,t) is odd (even) with respect to z for 
s £ T andte [0,T]. 

Remark 4.4. Note that if A(-, •, •) <G C°° (IxTx [0, T]) then the solution w{-,-, ■) of the problem P~2TJj) is 
a function in C°° (Kxfx [0,T]). This would be always the case whenever we apply Lemma\4.3\ 



The compatibility conditions must hold for every r in M., and so, in our derivation we refer to V := T t 
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4.2.1. Zero order term in e: For £ £ A/"(^f) using (|4.11[) wc obtain 

(4.26) = -fh"{z) + f(m(z)) for z £ [ - ^ 
while (gmi yields 

(4.27) = /(±1) for £ G fi \ «/V(e ). 

The above relations determine mo, more specifically mo{z) = m(z), where m solves the equation 

-m"(z) + f(m(z)) = 0, zel. 

4.2.2. First order term in e and determination of Vq via the Fredholm alternative. As explained at the 
beginning of this section, it is convenient for solving (|4.13[) to add a term eai(s, T)ifi'(z), s £ T and z£l, 
with ai(-,T) to be determined. This term will be subtracted to the second order. Recalling the definition 
of C, see (|4.2ip . and adding eai(s)rh' (z), we write (|4.13[) as 

(4.28) ^o{^ z , s) — f'(rh(z))(j)i(ez, s) — K(s)rfi'{z) + ea\(s)rh' (z) + G2fi{ez 1 s) = (£hi)(z, s), 
for s £ T, \z\ < §§. By (|4~14l) wc obtain 



(4.29) 



i(0 



/'(m(±oo)) 



eefi\AA(e ). 



We extend this definition of (f>i globally in 51. We then insert (|4.29|) into (|4.28[) obtaining for s £ T, 



(4.30) Mo(^,s) 



/'(™(*)) 
/'(±1) 



1 - 



/'(±1) 



G2,o(e z , s) — K(s)m'(z) + £ot\{s)m' {z) = {Lh\){z, s). 



Since the left hand side of (14.301) tends exponentially to as z — > ±00, if the solution of (|4.30[) exists (cf. 
Lemma |4~3|) . then it decays exponentially fast to 0. We can, therefore, extend (|4.30|) for any z in K. 
We prove the next result. 



Lemma 4.5. Set 



(4.31) 

where 
(4.32) 



V Q {0 \t,r)=Tr 



S<K(-)- 



J r K(QdS s \ 

|r| J 



+ 4 \ B ,o,o(-) ~ 



/ r B ,o,o(Od^ 



(0 



(m'(z)) 2 dz, 



and i?o.o,o(') *s defined in (|4.39[) lii/iife 7r is tte Dirichlet-Neuman operator given by (|6.4[) . T/ien i/iere 
exists a uniquely determined ai(-,r) £ C°°(r) and a unique solution h\(-,s) of (|4.30[) wii/i s £ T, swc/i £/iaf 
/ii(0, s) = and ft.i(-,s) £ L°°(R). Moreover, hi(-,s) for s £ T and its derivatives with respect to z decay 
exponentially fast to as z tends to ±00. Further hi(-,s) = hi(-,s) +eq £ (- 1 s) where hi(-,s) is an even 
function of z, /ii(0, s) = and h\(-, s) £ L°°(R). 



Proof. We start determining hi as solution of 



(4.33) no{0, 3 ) 



/'(±1) 



/'(±1) 



G 2 ,o(0, s) - J fC(s)m'(z) + ea^s)™' (z) = s). 
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Equation f|4.33|) differs from (|4.30[) only for terms of order e. Namely |G2,o(0, s) — G 2i o(ez,s)| < Ce\z\ and 
|/xo(0, s) — /j,q(sz, s)\ < Ge|z|. For any fixed sSf, the condition for the existence of hi requires that 



(4.34) 



Mo(0,s) 



1 - 



/'(±1) 



m'(z)dz + / G 2 , (0,s) 



/'(™(*)) 
/'(±1) 



m'{z)dz 



[K{s) -eai(s)] / (rh'(z)) dz for s G I\ 



Taking into account that f R f'(m(z))m'(z)dz = /(l) — /(— 1) = and J" m'(z)dz = 2 formula (|4.34[) can be 
written as 



(4.35) 2[/i (0, s) + G 2 , (0, s)] = [A(s) -eai(a)] / (m'(«))dz for seT. 

Recalling (|3.18[) we obtain 

2/*),o(0, s) = -2^0(0, s) - 2G 2 , (0, s) 



+ [if (s) - eai(s)] / (m'(z)) dz for s G I\ 



Mo,o,o(0 



y o (0) (r/)G(e,r,)d^+coW £ G fi, 



(4.36) 

Further, we set 
(4.37) 

(4.38) 

and 
(4.39) 

It is immediate to see that 
We first choose Vq ^ imposing for s G T, cfr (|4.36[) . the next identity 



Ao,o.o(0 = 2(y )r / G(Z, V )dS v - [ G(£,ri)G h0 (ri)dr) 
Jr Jn 



B ,o,o(s) = -2 [Ao,o,o(0, a) + G 2 , (0, s)] 

Mo,o(£) - Mo,o,o(0 - e, 
Mo(0 - £o,o,o(£) - £■ 



(4-40) 2/i ,o,o(0 = B ,o,o(0 + # (0 / ("*'(«)) dz CeT. 

Inserting (|4.37|) in (|4.40[) and integrating over T we obtain that 



(4.41) 
and therefore 



4 f dS e f V m (r 1 )G(Z, V )dS v + 2co(t)\r\ = f A(£)dS c [ (m'(z)fdz+ f B , ,o(£)d% 
Jr Jr Jr Jr Jr 



(4.42) 



K(QdS $ I {m'{z)Ydz + I S ,o,o(e)d^-4 / dS 6 I V {0) ffldS, 



We note that Co(i) is written in terms of the velocity field Vq which still needs to be determined. We insert 
Co(t), as in (|4.42l) . into (|4.40p to obtain the equation determining Vg ^ 

1 

4 

Here, Sr is the linear operator defined in (|6.5[) and 5 the quantity defined in (|4.32|) . Applying the Dirichlet- 
Neumann operator, see (|6.6p . we arrive at f|4.31|) . The above determines first Vq " 1 and then co(i), see (|4.42[) . 
Now since Vq and Co(t) are chosen, we may then simply choose &i(s) so that (|4.34p is satisfied. Then for 



S r V^( V ,T) = S 



T3 / r Bo,o,o(C)d^ 
£>o,o,(U5j 



iri 



rier. 



J r K(QdS, 

k(v) — 
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any sgT, the existence of a unique hi satisfying /ii(0, s) = is assured; h\ is exponentially decaying to zero 
as \z\ — > oo. Since the left-hand side of (|4.33[) is even, then the solution hi(-, s) is even as a function of z. 
Define eq e (-, s) = hi(-, s) — hi(-, s) and subtract (|4.33[) to (|4.30|) . We have that q e (-, s) satisfies 



(4.43) 



1 - 



/'(™(*)) 
/'(±1) 
(Cq £ )(z,s), 



[li (£z,s) - /Lt o (0, s) + G 2 fi{ez,s) - G 2j0 (0,s)] + [ai(s) - ai{s))m'{z) 



where ai(-) must still be determined. We determine ai(-) so that the following solvability condition for 
q e (-,s) holds: 



1 



— / dzm'(z) 



1 



/'(±1) 

+ [ai(s) - &i(s)] I dz(m'(z)) 2 = 



\fM)(ez, s) - /i (0, s) + G 2 fi{ez 1 s) - G 2 ,o(0, s)] 



(4.44) 



By Lemma T4.3l we have that for all s £ T, q e (0, s) = 0, q e (Q, s) E L c 
Remark 4.6. Let us denote by 



□ 



A*o,o,o — Mo, o,o + Ao,o : o> 

the quantities defined in (j4.37[) and (|4.38[) . Taking into account (|2 . 19[) and the definition of (Vo)r = 
JprTJ fnGifiiv)^; see (|3.17}> . we /lave i/iaf /i^ ,o * s ^ e solution of 

(4.45) A M (0 = -Gi,o(0 for £ e fi \ T, 

(4.46) /x(£) =2S , ir(C)-G 2 ,o(0 on I\ d v (i = on 90. 
One deduces (|4.46j) &?/ (|4.40[) , taking into account the definition of S and (|4.39|) . 



4.2.3. Second order term in e. From (|4.16|) we have that 
1 



(4.47) 



>(0 



/'(±1) 



M0 - / 2 (m(±cx))^i(e)) + eA& + G a ,i(0] , £ e n\jV(e ). 



As was done before, we extend the validity of (|4.47[) globally in Q. We insert (|4.47|) into (|4.15[) . Further, we 
add and subtract to the next order ea 2 (s)fh' (z) to obtain 



(4.48) 

where we set 



fii(ez, s) 



/'(±1) 



A 2 (z, s) + ea 2 (s)m'(z) = (Ch 2 )(z, s), 

f'(rh(z)) 



(4.49) 



A 2 (z,s) = [G 2 ,i(ez,s) + sA<f> 1 (ez,s) - f 2 (±l,(j>i)(ez,s)] 
- K 2 (s)zfh'(z) - K(s)h' 1 (z,s) + a 1 (s)m'(z). 
All the quantities in (|4.49[) have been already determined. Furthermore, 



/'(±1) 



(4.50) 



lim A 2 (z, s) = s G T, 

|#|— >-oc 



the convergence being exponentially fast. As done before, we extend 1)4.480 in R. 
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Lemma 4.7. Set 

(4.51) 



y 1 (0) (c,r) = r r 



(0 



where Bi(s) is defined in (|4.54p and Tr is the Dirichlet-Neumann operator. Then there exist uniquely 
determined a 2 (-, T) G C°°(T) and h 2 {-, s) G L°°(R) with h 2 (0, s) = for s G T, solutions of (|4.48j) . Moreover 
h 2 (-,s) and its derivatives with respect to z decay exponentially to 0, as z tends to ±00. 

Proof. The solvability condition, see (I4.24[) . is satisfied provided that for s G T and t G [0, T] the next relation 
holds true 



(4.52) 



Hx{ez, s) 



/'(I) 



m'(z)dz= / A 2 (z, s)m / (z)dz — ea 2 (s) / (m'(z)) dz 



where fj,i is defined in (|3.34j) . The term fix of fi\ has been already completely determined. As in the previous 
case, still to be determined are the constant ci(t), the velocity V{ and a 2 (s). First, we write (|4.52j) as 



(4.53) 
where 
(4.54) 

Let now 
(4.55) 



f , ,\ f'(m(z)) 



Bi(s) 



Mi,o,o(£) 



A 2 (z, s) - fii(ez, s) 



m'(z)dz = Bi(s) — ea 2 (s) j (m'(z)) dz, 
f'(m(z)y 



1 - 



/'(I) 



m'(z)dz. 



Since /ii,o,o(£) — Mi,o(£) — e : we choose V\ by imposing 

/'(m(z)) 



(4.56) 

and obtain 
(4.57) 



A*i,o,o(0,s) 



/'(I) 



m (z)dz — Bi(s), 



Mi,o,o(0, s) = 2 B i( s ) 



s g r. 



Inserting (|4.55|) in (|4.57p and integrating over T we arrive at 



(4.58) 



\( B 1 {r,)&S n -2 f AS J v} 0) (r,)G({, V )dS n 



Observe that L (s)ds = and let Sr be the linear operator defined in (|6.5j) . Then obviously, (|4.57[) can 
be written as 



SrV 1 {0 \0 = \B 1 (0-^r l J r B 1 (s)ds feT, 



and thus, applying the Dirichlet-Neumann operator (see (|6.6|) ) we obtain (|4.51|) . This determines the (con- 
stant in £) ci(t). Now, since Vi (0) andci(t) are determined we may choose a 2 (s) so that (|4.53|) is satisfied. □ 

Remark 4.8. Notice that /ii,o,o solves 

A/i 1A o = o for £Gft\r, 

(4.59) 1 
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4.2.4. n—th order term in e, 3 < n < N : As previously, we determine the function <p n for £ G fi \ A/"(eo) 
from ()4.18j) . Then, we extend the validity in f2 obtaining 

(4.60) MO = J^[^n-l(0+^n-l(0-fn(±l,4>l,<t>2,--- , <^»n-l)(0 + G 2 ,n-1 (0] £ G 

We then insert (|4.60[) into (|4.17|) . we add and subtract to the next order the quantity ea n (s)fh'(z), to the 
left hand side of (|4.17[) and we obtain 

f(m(z)Y 



(4.61) 
where we set 



/'(I) 



A n (z, s) + ea n (s)m'(z) = (Ch n )(z, s), 



n-l 



A n (z, s) = -a n -i(z,s)m' - [a„_j(z, s)/i-(z, s)] - b n ^ i (z,s)—h l (z,s 



; = 1 



n — 3 



c n -i(z,s) — hi(z,s) 
as 



eA(/) n _i (ez,s) 



1 - 



/'(i) 



1 f'(rn(z)) 
/'(I) 



/'(±1) 

It is easy to verify that for all s G T 



/n(±l,</>i,02, • • • I ^i)(K,s)-/n(mo I ™i 1 "i2r" , m n ^ 1 )(ez, s) . 



lim A n (z, s) = 0, 



where the convergence is exponentially fast. Namely there is no problem for those terms involving to', /ij(-, s) 
and their derivatives because of the exponential convergence to zero of all these terms for all s G T. For 
the remaining terms recall that lim /'(to(z)) = /'(±1), = hi + 4>i with /ii(z, s) — > as \z\ — > oo for all 

|z|— >oo 

s G r, all limits being exponentially fast. Then one obtains immediately 

f'(fn(z)y 
/'(±1) . 

exponentially fast also. We extend (|4.61|) to hold on all of R and regard it as an equation for h n (-,s) for 

s G r. 



lim 

\z\— >oo 



1 - 



/ n (±l,01,02, ' ' ' ,4>n-l){£Z, s) = 0, 



Lemma 4.9. For any positive integer n, n < N , set 



(4.63) 



V^T)=T r 



1 



B„_!(s)ds 



/or e G T, 



where B„_i(s) is defined in (|4.66[) . T7ien i/iere exist uniquely determined a n (-,T) G C°°(r) and h n (-,s) G 
i°°(R) /or s G r wii/i /i„(0, s) = 0, solutions of (|4.61[) . Moreover, h n {-,s) for all s G T, and its derivatives 
with respect to z decay exponentially to as z ±oo. 

Proof. The solvability condition is satisfied provided that 

/'(m(z)) 



(4.64) 



1 



/'(I) 



m' {z)dz — \ A n (z,s)m'(z)dz — ea n (s) / (to'(z)) dz 



Since in view of Q3.38P /j„-i = fJ-n-ifl + P>n-i with /2„_i being already determined to satisfy (|4.64[) . we 
require that 



(4.65) 
where 
(4.66) 



fi n -i fi (ez,s) 



1 



/'(I) 

A n (z,s) - s) 



m'(z)dz = B n -i(s) - ea n (s) / (m'(z)) dz, 
/'(to(z))" 



/'(I) 



m'(z)dz. 
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^-1,0,0(0 = 2 / V n - 1 ( V )G^,r ] )dS v +c n ^ 1 {t). 



We set 
(4.67) 

Since /i„_i j o(^, t) — (J, n -i,Ofl(£-,t) — £ then we may determine V n -\ by imposing 

f{m{z))- 



Mn-l,0,o(0,s) 



1 - 



/'(I) 



m'{z)dz = B n _i(s), 



//„-i,o,o(0,s) = ^-B n _i(s) 



s g r. 



(4.70) 



(4.68) 

obtaining thus 
(4.69) 

Inserting (I4.67P in (|4.69f) and integrating over T we get 

1 

IT 

We insert (|477D|) into (|437|) . so, ([41)5)1 gives 

and (|4.63[) follows. This determines V^°\ first and then c n -i(t). Therefore, we may define a„ in order to 
satisfy ([OS)) . □ 

4.3. Proof of Theorem 14. li To complete the proof of Theorem 14.11 we need to estimate the remainder 
term (cf. (|4.19l) ) given by 



(4.71) 



£ i? 2 (C, i, e) = eR 2 (Z, t, e) - e ^a N (s(0,t)m'(^^-). 



Since (|4.20l) holds true, we obtain that 
(4.72) 



sup sup \Ra(Z,t,e)\ < Ce N . 
fen te[o,T] 



So, Theorem 14. II is proved. 



□ 



5. Proof of Theorem 12.51 
The proof of Theorem 12.51 is an immediate consequence of the following result. 



Theorem 5.1. Take N > 1 and G\ and Gi as in Assumption Al. There exist vector fields Vj, j = 
0, • • • , (N — 1) and functions rrij, j = 0, • • • ,N from Ai to C°°(0) as in the Ansatz such that the 
following holds. For any Tq G Ai, choose ko > fc(ro), set £0 = ^ and let T be the lifetime of the solution of 
(|2.7[) in M., according to (|2.3[) . Then for all t < T and for all e G (0, £0] we can construct (fh^ N \ ft^ N ~ 1 ') G 
C°°(nx[0,T}) where mW is an^ modification ofm^ N \ i.e. sup (et)enx[0T] |m( w )(f,t)-mW^i)| < Ce N 

and ft^-V is an e"- 1 modification of n (N ~ 1} , i-e- sup (f)i)enx[0) ' T] |/U (Jvll) (£,*) - A* (JV-1) (6*)I < Ce^ 1 
satisfying 



N-l 



(5.1) 



JV-1 



t ) = - £ AmW(e,t) + -/(m W (£,*)) - ]T ^G 2 , ■,-(£) + fl^K, i, e) *n fi x (0,T), 



3=0 
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where 



sup sup 

£entg[o,T] 



< Ce 



N-l 



Further, fL^ N ~ x ' (-,t) and rh^ N \-,t), fort <G [0,T], satisfy Neumann homogeneous boundary conditions on 
the boundary of CI. In addition 

,F 



(5.2) 

where [1q is the solution of 
(5.3) 

subject to the boundary conditions 



sup sup|M (JV - 1) (£,t)-/xSo,o(£ ) *)l <Ce, 
te[o,T] £efi 



A/i(0 = -Gi >0 (0 C^\r 



(5.4) 
and 

(5.5) 
(5.6) 

Proo/. Set 

(5.7) 

where 



fi(C) = 2SJr(0-Ga,o(0 on if" , 5^ = on 90 



sup sup 

te[o,T] CeA A(e ,r< N) ) 



mW(e,t)-m 



'erf; 



< Ce 



sup sup 

*eto.3l fgn vv(6o),rW) 



m W (e<)Tl 



mW^^mWf^)- / R 1 (r,e)dr, 

Jo 

1 
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i*i(£,*,e)d£, 



(5.9) 



and i, e) is the remainder in Theorem 13. 11 defined in (|3.14[) and estimated in p.3[) and (|3.4[) . Let us 

denote by 

(5.8) ^ N -V(U)=V {N - 1] (U)+v(U), 
where v(£,t) solves 

Av(^t)=Rx(i,t,e)-R 1 (t,e) for ^£fi, 
9„v = on dfl, 
with the further requirement 

f v(£,t)d£ = o te[o,T]. 

Jn 

Since \Ri(£ t t,e)\ < C(T)e Ar " 1 we have that \v(£,t)\ < Ce N -\ The functions # and /Z^" 1 ) satisfy (1531) . 
Namely the first equation of (|5.1[) is satisfied by Theorem 13.11 and by construction, see (|5.7[) and (|5.9p . The 
second equation is obtained from Theorem 14.11 adding and subtracting terms to obtain p,^ N ~^ and m^ N '. 
We obtain 



where 



£(AT-1) = ^(Ar-1) + „ = _ £Aj ^(iV) + i /( aW) + 

e 

JjW=JjW(f,t, e ) = i /(mW+| t E 1 (T )£ )drj-/(mW) 



i?2 + U, 



and i? 2 is the remainder in TheoremED see (j4TT|) . Since Pi = ©(e^), P 2 = ©(e^), w = ©(e^ 1 ) and 



(5.10) 



/(mW+ / Pi(r,e)dr I - f{rh (N) ) 



< 



P 1 (r, £ )dr-0(e JV - i ), 
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then the second equation of (|5~T|) is satisfied as well. In Remark H~6] it is explained that if ^ 

0,0 — Mo, 0,0 + 

Ao,o,Oj where /J<o,o,o and Ao.0,0 arc the quantities defined in (|4.37l) and (|4.38[) then it verifies (|5.3[) and (|5.4[) . 
The relation (|5.2p is then immediate from their definition and Theorem 13. II The (|5.5[) and (|5.6|) are satisfied 
by construction of the mS N \ Theorem 12.51 is then proved. □ 



6. Appendix 

A.l: The Dirichlet Neumann operator. We recall in this section the main properties of the Dirichlet- 
Ncumann operator which we have been using through the paper. Some of these results were already presented 
in the Appendix of [TT|. Let G(£,,rj) be the Green's function in f2 defined in (|2.19p and verifying (|2.20[) . To 
define the Dirichlet Neumann operator we consider the following single layer potentials. Given a smooth 
function h defined on F G Ai , consider the single layer potential 

where dS v denotes the arc length measure along T. The function 0^ satisfies a Neumann boundary condition 
on dCl, and also the equation 

a<M£) - h(0 - p J %)dS„ £ g r. 

The curve T separates O in two subsets. We denote by fip the interior of T and by Hp" its exterior. We 
denote by n the unit outer normal to i7p . There is a discontinuity in the normal derivatives of <ph across T, 
so, we have that 

(6.1) no = \ \bm t (0, 

where the right-hand side denotes the jump of the normal derivatives at £ G T, i.e. 

[dvMv (0 = (Wn+ (0 - (d.Mn- (0- 
This is a well known result from potential theory |21j . For £ away from F, 

a^(0 = -pr J r H v )ds n . 

Thus, the single layer potential is harmonic away from T if and only if J r h(£)dS£ = 0. Otherwise, it is 
subharmonic or superharmonic, according to whether — J r h(S,)dS^ is positive or negative. Every continuous 
function <fi harmonic away from T, satisfying the Neumann boundary condition and the following relation 

(6.2) f ^ = 0, 

is the single layer potential of a uniquely determined function h defined on T and satisfying 

(6.3) J h(£)dSt = 0. 

Indeed, if <fih is such a single layer potential, then from (|2.20[) , we get that J Q </>h(£)d£ = 0. 

On the other hand, let <fi be any continuous function that is harmonic on fip and f2p~ , and which satisfies 
(|6.2p . Let us define h in T by 

MO = \ [Mr (0> 

and refer to this as the Neumann data for <f>. By the- divergence theorem we obtain 

2 f /i(£)dS £ = f S^ + d5 £ - [ d v </>-dSt = - [ A<£df = 0, 
Jr Jr Jr Jn\r 

where ± denotes the restriction of <f> in f^. Hence, h satisfies (|6.3[) . 
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Notice that <f> — fa satisfies the Neumann boundary conditions and 

[W - &)] r (0 = o. 

This means that (j) — fa is a constant. Since the integral is zero then <f> = fa. This proves the one to one 
correspondence between single layer potentials of functions h satisfying (|6.3p . and continuous functions <f> 
that are harmonic on Sip and Sip", and satisfy (|6.2[) . 

Next, given a continuous function <f> that is harmonic on Sip and Sip , whether or not (|6.2[) is satisfied, we 
define the function g on r by g := <j)\r- We naturally refer to g as the Dirichlet data for cf>. The Neumann 
data is [d v 4>] T . The Dirichlet-Neuman operator 7r is defined by 



1 



(6.4) Tvg = - 



Jr • 



where <f> is the continuous function that is harmonic in Sip and Sip , with </>|r = g- 

A simple argument shows that 7r is a positive Hcrmitian operator. Indeed, let ip be continuous on SI, 
and harmonic on Sip and Sip , with ?/>|r = h. Then it follows that 

2 J hTvgds = J ^[d v <t>] T ds 

Op J Op 

VV> ■ V0d£. 

n 

Taking ft, = 1, so that ip = 1, we further see that the range of 7r is orthogonal to the constants. We let 7r 
denote the Friedrichs extension of 7r- It is easy to see, and well known, that the form domain of 7r is the 
Sobolev space i? 1 / 2 (r), and that the kernel consists exactly of the constants. There is an explicit formula 
for the inverse of 7r restricted to the orthogonal complement of the constants; we denote this by Sr- Indeed, 
let v be any function on T with J r v(s)ds = 0. Since the single layer potential <p v for v has Neumann data v, 
all we need to do is to subtract a constant to make this function orthogonal to the constants on T, instead 
of being orthogonal to the constants on SI. Therefore, the inverse Sr is given by 

(6.5) S r v(0= I G(tvHv)dS v ~^- [ [ G(tvHv)dS v dS e ^T. 



It is easily checked that the inverse operator is self adjoint on the orthogonal complement of the constants. 
Now let h be an arbitrary smooth function on T satisfying J T h(s)ds = 0, and consider the single layer 
potential 

0(e) = ^G(e,r?)%)d^ £6 0. 

In general, the Dirichlet data for <fr do not integrate to zero on T and hence, are not directly related to the 
Neumann data through the Dirichlet-Neumann operator. However, we can correct this by subtracting a 
constant and defining the function 

0(0 = ^(0-4 / m*s n . 



Then obviously we have 

0|r 



= S r h, 



(6.6) 

h = Tr(p 



We can now express the vector field V driving the Mullins-Sekerka flow as 

1 

W\ 



(6.7) V = Tr ( K - / K(s)ds 



r 
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We close by establishing notation for the two harmonic extension operators that will arise throughout 
what follows: 

The Neumann harmonic extension operator £r,N is defined by 

(6.8) (£r, N v)(0 = J G^v^dSr,-^! J G^v^dS^ £ e 0, 
where v is a function on F satisfying 

y«(f)dS c = 0. 

Notice that (£r,NV ) (£) is the unique function that is continuous on f2, harmonic on Q\T satisfying Neumann 
boundary conditions on dfl, with Neumann data v, and with zero integral over F. 

The Dirichlet harmonic extension operator £t,d is defined by setting £r. to be the harmonic function 
(j> on f2\r with Neumann boundary conditions on dT, and 0|p = g. Here, there is no restriction on the integral 
of g over T. Naturally the Dirichlet extension can be expressed in terms of the Neumann extension and the 
Dirichlet-Neumann operator. Relations (|6.5[) and (|6.8[) give that 

(6.9) £ rD5 ( € ) = £ rN ^rr^-py5Wd^))(0 + py.9Wd^ £ e 0. 



A. 2: The expansion in e of the Laplacian in local coordinates. Let f(z,s), with z = |, be a C 2 
function from R x T to R. Then, in the two dimensional case, we have that 



e 2 A/(z, s) = — — ((1 - K(s)sz)f z ) z + s 



1-K(a)ez\" w \l-K(s)ez 
(6.10) k V j 

l « i~si q £k(s)z 



1 - X(s)ez (1 - K(s)ez) 2 J (1 - if(s)ez) 3 ' 
Recalling that for |x| < 1 it holds that 



^ oo oo 1 1 °° 

v ; n=0 v ' n=0 v y n=0 

we may rewrite (|6.10[) as follows 

oo 

(6.11) e 2 Af = f zz + J2 Wn+i(z, s)f z + b n+1 (z, s)f ss + c n+1 (z, s)f s } , 



n=0 



where 



3.12) 



a n+1 (z,s) = -K n+1 (s)z n , 
b n+1 (z,s) = nK n -\s)z n -\ 

c n+ i(z,s) = \n{n- l)z n - 1 K n - 2 (s)^K(s). 
2 as 

A. 3: proof of Theorem 12.71 The proof goes very much as in [IJ Theorem 2.1]. Hence we outline only 
those points where the presence of the G\ and Gi makes a difference. First of all the constructed functions 
m( N ^(t), t £ [0,T] satisfy the requirements needed to apply the spectral estimates proven by [2] and |12j . 
Namely, sec (|2.8|) . (|2.10[) and Lemma l4~5l we can always write for £ € f2 and £ [0,T] 
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where mo is given in (|2.9j) . h±(-,s) is the function determined in Lemma 14.51 which is even as function of 
2 € 1 for any s g T { t N) , equal to in n\Af(e ) and when d(£,,T { t N) ) = 0. We denote by e 2 q £ (£,,T { t N) ) the 
remaining functions in the expansion of rrS N ^ which are equal to zero in O \ N{eq) and by £0 e (£, ) the 
corrections to ±1 in Q \J\f(eo). Recall m^ N \-,t) arc C°°(fi) for any t G [0, T]. We immediate have, since 
m(-) is odd while hi(-,s) is even, 



hi{z, s){m'{z)) 2 f"(fh(z))dz = 6 / tn(z, a)(m' \z)) 2 fh{z)dz = 0, Vs e T\ 



(N) 



This is one of the requirement needed to apply the spectral estimates. The remaining requirements are 
immediately satisfied by the smoothness of m^ N '(-,i) and by the fact that the $j, j = 1,...,N in the 
expansion of m^ N > satisfy a global Lipschitz bound independent on e. 
Then one proceeds as in [I]. Write (|T8)) as the the following : 



(6.13) 



dtm 6 = A/i e + Gi in Qx, 

u E = — eAm £ + -f(m £ ) — G2 in Qt, 

£ 

m £ (£,0) =m (£) £e O, 
d„m e = d v /\m £ —Q on <9S1, 



and (|2.12j) the the following: 

in Vtr, 



N-l 



3=0 



N-l 



3.14) 



AN) 



-eAmW + - /(mW) - ]T ^G 2j + flW in T 



cW^ = ^Am^' = on 9fi. 

Define R := m £ — rrS N \ Then integrating i? in space, by (|6.13[) and f|6 . 14[) and the fact that 0) = 0, we 
obtain 



m e_ m (N) U£ 



(6.15) 



d^J d s (m £ -m (JV) )ds 
/ d£ / (A// + d - A/x (iV ) - V e>GjJ)da 



/ A 



r/.s 



(AT) 



(AT) 



d£ = 



N 







fdsf 


Gi.N 


10 Jn 





d^ds 



since d u [fi £ — fi^ N '] — on dfl, and (|2.2ip holds. Note that we need (and used for the above) d v G2 = 
d u R( N ^ = on <9f2 in order to have d v [fi £ — /■« ] = on dQ. In addition, we need the residual R^ 
satisfying a Neumann condition i.e. 

<9„i? (Ar) = on dn 

which is true by construction, cf. system (|5.1[) and the b.c. on system solutions in the statement of Theorem 
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Hence, since by (|6.15[) . for any t € [0,T] J n i?(£,t)d£ = then there exists unique ip(£,,t) such that 

- AiP{-,t) = R(-,t) in il, 
d v ij){-,t) = on dfl, 
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